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OPTICAL POLARIZATION AND SPECTRAL VARIABILITY IN THE M87 JET 
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Raymond C. Simons 1 , Eric Clausen-Brown 7 , C. C. Cheung 8 , Lukasz Stawarz 9,10 , Markos Georganopoulos 11 , 

William B. Sparks 12 , John A. Biretta 12 

ABSTRACT 

During the last decade, M87’s jet has been the site of an extraordinary variability event, with one 
knot (HST-1) increasing by over a factor 100 in brightness. Variability was also seen on timescales 
of months in the nuclear flux. Here we discuss the optical-UV polarization and spectral variability of 
these components, which show vastly different behavior. HST-1 shows a highly significant correlation 
between flux and polarization, with P increasing from ~ 20% at minimum to > 40% at maximum, 
while the orientation of its electric vector stayed constant. HST-l’s optical-UV spectrum is very hard 
(auv-o ~ 0-5, F v oc and displays “hard lags” during epochs 2004.9-2005.5, including the peak 
of the flare, with soft lags at later epochs. We interpret the behavior of HST-1 as enhanced particle 
acceleration in a shock, with cooling from both particle aging and the relaxation of the compression. 

We set 2a upper limits of 0.55 parsecs and 1.02c on the size and advance speed of the flaring region. 

The slight deviation of the electric vector orientation from the jet PA, makes it likely that on smaller 
scales the flaring region has either a double or twisted structure. By contrast, the nucleus displays 
much more rapid variability, with a highly variable electric vector orientation and ’looping’ in the 
(I, P ) plane. The nucleus has a much steeper spectrum ( auv-o ~ 1-5) but does not show UV-optical 
spectral variability. Its behavior can be interpreted as either a helical distortion to a steady jet or a 
shock propagating through a helical jet. 

Subject headings: galaxies: individual (M87) - galaxies: active - galaxies: jets; nuclei 


1. INTRODUCTION 

M87’s jet was one of the first manifestations observed 
of the active galactic nucleus (AGN) phenomenon (Cur- 
tis 1918), and has been the target of myriad observations 
due to its brightness and also proximity ( d = 16 Mpc, 
Tonry 1991). During the last decade, M87’s jet has been 
the site of an extraordinary variability event, with knot 
HST-1 increasing in optical/UV brightness by a factor of 
more than 100 between 2000 and its peak in 2005. The 
flare in knot HST-1 has been the target of several moni- 
toring efforts using the Hubble Space Telescope (hereafter 
II ST), Chandra X-ray Observatory , VLA and other tele- 
scopes. Previous papers from this project include Paper 
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I reporting our first results (Harris et al. 2003), Paper II 
which focused on the HST data (Perlman et al. 2003), 
Paper III which was mainly on the X-ray lightcurve of 
HST-1 which delineated the massive 2005 flare (Harris 
et al. 2006), Paper IV which focussed on the VLB A 
results, showing superluminal proper motions in HST-1 
(Cheung et al. 2007), and Paper V (Harris et al. 2009) 
which focused on a more detailed analysis of the variabil- 
ity timescales of HST-1 and the nucleus. Madrid (2009) 
has also added a complete analysis of the UV light-curve 
of HST-1 and the nucleus between 2000-2006. 

In this paper (VI of the series), we discuss two addi- 
tional aspects of the monitoring campaign, namely the 
evolution of the polarization and spectral index in the 
optical-UV. We concentrate on the nucleus and HST-1, 
as they are the main variable components in the jet. A 
future paper will combine these observations to produce 
a new polarization map of the entire jet and discuss any 
changes over the decade between the data of Perlman et 
al. (1999) and this paper. Section 2 provides a detailed 
explanation of the observations and data reduction pro- 
cess. Then, in Section 3, we will discuss the techniques 
used in reducing the data. In section 4, we will discuss 
results. Finally, in Section 5, we discuss implications for 
jet variability models and close in Section 6 with a sum- 
mary. 

2. OBSERVATIONS 

As has been well documented (e.g., Paper I, Waters 
& Zepf 2005), the flare of knot HST-1 began sometime 
during 2000. While the jet of M87 was a regular target 
for HST and Chandra almost from the start, intensive 
monitoring by these telescopes began in 2002 (Papers 
I, II, III, Madrid 2009). Here we review these observa- 
tions, concentrating on the optical polarimetric part of 
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the campaign. 

Optical polarimetry was obtained on a somewhat dif- 
ferent schedule than the UV imaging discussed in Madrid 
(2009) Of the eighteen observations obtained between 
2002 December and 2007 November, fourteen occurred 
on the same schedule as the imaging observations during 
November 2004- December 2006. The other four obser- 
vations were at roughly yearly intervals before and after 
this period. Polarimetry was done in two bands, F330W 
and F606W, with F606W observations being done much 
more often. Table 1 details the scheduling of these obser- 
vations, along with other information about the nucleus 
and HST-1 which will be described later. 

The High-Resolution Channel (HRC) of the Advanced 
Camera for Surveys (ACS) was used to obtain the po- 
larimetry data for 17 of the 18 epochs. The ACS 
HRC is a single-chip CCD camera, with a plate scale 
of 0.028x0.025 arcsec/pixel, corresponding to a field of 
view of about 28 x 25 arcseconds and yielding diffrac- 
tion limited resolution of «0”.06 for the F606W obser- 
vations, and a Nyquist-limited resolution of % 0” .05 for 
the F330W observations (see Maybhate et al. 2010). On 
ACS the polarizing set is comprised of either the POLOV, 
POL60V and POL120V filters (used for F606W) or the 
POLOUV, POL60UV, and POL120UV filters (used for 
F330W). As the three polarization filters can be selected 
individually Maybhate et al. (2010), no change in posi- 
tion or chip was necessary between them. All observa- 
tions were CR-SPLIT to mitigate the effects of cosmic 
rays, but dithering was typically not done. 

For the final epoch, which occurred after ACS went 
offline due to an electrical short, the Wide Field Plan- 
etary Camera 2 (WFPC2) was used for polarimetry. 
On WFPC2 the polarizing filter set is comprised of the 
POLQ quad, which has filters at angles of 0, 45 and 90 
degreees. The WFPC2 is a chevron-shaped, four- CCD 
camera, with three wide-field chips (WF2, WF3, and 
WF4) and a fine-resolution one (PCI). Its setup is less 
flexible for polarimetry, as the POLQ quad can only be 
rotated through 51 degrees (Biretta & McMaster 1997). 
For this reason, we used the three WF chips. These chips 
have a plate scale of 0.09965 arcsec/pixel, leading to a 
final resolution ~ 0.2” set by the Nyquist theorem. 

3. DATA REDUCTION 

We obtained the data for these observations from 
the HST archive. All data were recalibrated with up- 
dated flat field files and image distortion correction 
(IDC) tables, obtained from from the STScI Calibra- 
tion Database System. Standard techniques were used 
to recalibrate the data. These methods are described in 
detail in the instrument handbooks for ACS and WFPC2 
(Maybhate et al. 2010; Mobasher et al. 2002). CTE 
corrections were computed using the data found in the 
ACS Instrument Handbook (Maybhate et al. 2010). 

After the data were retrieved, MULTIDRIZZLE was 
used to combine and cosmic-ray reject the images. The 
alignment of images was refined, assuming the positions 
of the core and HST-1 as canonical and using TWEAK- 
SHIFTS in IRAF to shift the images to a common frame 
of reference and correct for geometric distortions using 
the models in the IDC. This was done because experience 
with the polarizing filters on both ACS and WFPC-2 has 
shown that there can be small irregularities in the filters 


that can cause sources that, are far away from one’s re- 
gion of interest to shift apparently as compared to other 
sources in the field. The procedure was also checked by 
using only the nucleus as canonical, thus allowing for 
the possibility of motions in HST-1. Unfortunately with 
the short exposure times (typically only 600s per polar- 
izer) there were few or no globular clusters that could 
be used for all frames, especially as the chip orientation 
changed from epoch to epoch. Extensive testing gives us 
confidence, however, that our procedure successfully and 
repeatably aligns the images to ±0.2 pixels. Following 
MULTIDRIZZLE, the orientation was set so that the y- 
axis corresponds to north. The final result is a cosmic 
ray rejected and geometrically corrected image for each 
polarizer at each epoch (Fruchter and Sosey 2009). 

For the F606W data, it was necessary to subtract 
galaxy emission before performing photometry and po- 
larimetry. We first created a composite image of the 
galaxy by drizzle-combining all epochs together, to im- 
prove the S/N on the host galaxy. We then modeled the 
galaxy emission using the ELLIPSE and BMODEL tasks 
in IRAF. After the model image was computed, it was 
split into 3 models to correspond to each polarizer, and 
subtracted, using the IRAF command IMCALC, from 
the corresponding drizzled image for each polarizer. 

3.1. Polarization Images 

Next, the drizzled images from each polarizer were used 
to create images for Stokes I, Q, and U, along with their 
errors. For the ACS data, we followed the procedure in 
the ACS data handbook (Maybhate et al. 2010). For 
the WFPC2 data, the Stokes images were computed by 
using the WFPC2 Polarization Calibrator tool 13 . This 
produces the coefficients needed to cbmpute the Stokes 
images by using Mueller matrices that describe the pick- 
off mirror, the polarizer filter, and the various rotations 
between the optical elements and the reference frames. 
This tool is accurate to ssl-2% (Biretta & McMaster 
1997). Both of these procedures yield Stokes U , Q and 
I images that are combined in a standard way to pro- 
duce emission weighted fractional polarization (defined 
as P = ( Q 2 + U 2 ) 1 / 2 //) and electric vector position an- 
gle (defined as EVPA = 1/2 x tan -1 (C7/Q)) images. 

After the Stokes images are found, we accounted for 
the well-known Rician bias in P (Serkowski 1962) using 
a Python code adapted from the STECF IRAF pack- 
age (Hook et al. 2000). This code debiases the P im- 
age following Wardle &: Kronberg (1974), and calcu- 
lates the error in polarization PA, accounting for the 
non-Gaussian nature of its distribution (see Naghizadeh- 
Khouei & Clarke 1993). In performing this calcula- 
tion, pixels with signal to noise (S/N) <0.1 were ex- 
cluded outright, and since the debiasing is done with a 
“most-probable value” estimator, pixels where the most- 
probable value of P was negative, or above the Stokes I 
value ( i.e . P > 100%) were blanked. This code was first 
used in Perlman et al. (2006). 

3.2. Aperture Photometry and Polarimetry 

To obtain fluxes in the Stokes parameters, we per- 
formed aperture photometry. For the ACS data, we used 

13 http://www.stsci.edu/hst/wfpc2/software/wfpc2_pol.calib.html 
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TABLE 1 

Photometry of M87 Components 


Date (No.) 

F606W 

Core Fluxes (pJy) 
F330W F250W 

F220W 

F606W 

HST-1 Fluxes ( M Jy) 
F330W F250W 

F220W 

Dec 07 2002 

(1) 

671 ±7 

305 ±2 


146 ± 11 

232 ±2 

231 a ± 2 


141 b ±6 

Dec 10 2002 

(2) 

630 ±6 

305 ±2 


146 ± 11 

237 ±2 

231 ±2 


141 ±6 

Nov 29 2003 

<3) 

478 ±5 



137 ±11 

427 ±4 



244 ±7 

Nov 28 2004 

4) 

1066 ± 11 

475 ±2 

306 ± 14 

226 ± 14 

970 ± 10 

882 ±4 

674 ± 12 

631 ± 12 

Dec 26 2004 

(5) 

1306 ± 13 


363 ± 10 


1113 ± 11 


719 ±10 


Feb 09 2005 

(6) 

891 ±9 


280 ±9 


1224 ± 12 


738 ± 10 


Mar 27 2005 

(7) 

1037 ±10 


328 ± 10 


1404 ± 14 


903 ± 11 


May 09 2005 

(8) 

932 ±9 

446 ±3 

274 ±9 

217 ±9 

1333 ± 13 

1209 ± 5 

878 ± 14 

853 ± 14 

Jun 22 2005 

(9) 

839 ±8 


273 ±9 


1150 ±12 


877 ±11 


Aug 01 2005 

(10) 

639 ±6 


192 ±7 


1117 ± 11 


664 ±10 


Nov 29 2005 

11) 

735 ± 7 

349 ±2 

217 ±8 

170 ±8 

1019 ± 10 

872 ±6 

604 ± 12 

616 ± 12 

Dec 26 2005 

(12) 

756 ±8 


234 ±8 


987 ± 10 


590 ±9 


Feb 08 2006 

13) 

631 ±6 


201 ±8 

160 ±8 

771 ±8 


482 ±8 

469 ±8 

Mar 30 2006 

(14) 

780 ±8 


232 ±8 


656 ±7 


404 ±7 


May 23 2006 

(15) 

862 ±9 

372 ±2 

193 ±7 

156.4 ± 7 

592 ±6 

488 ±4 

361 ±7 

357 ±7 

Nov 28 2006 

16) 

1370 ± 14 

636 ±3 

323 ±9 


862 ±9 

720 ±4 

486 ± 10 


Dec 30 2006 

(17) 

1006 ± 10 


276 ±9 


682 ±7 


414 ±7 


Nov 25 2007 (18) 

1292 ± 17 




372 ±4 





a Observations taken Dec. 10, 2002. 
b Observations taken Nov. 30, 2002. 


TABLE 2 

POLARIMETRY AND SPECTRAL INFORMATION 


No. 

F606W 

P(%) 

F606W 

EVPA(°) 

Core 

F330W 

P(%) 

F330W 

EVPA(°) 

<*o-uv 

F606W 

P(%) 

F606W 

EVPA(°) 

HST-1 

F330W 

P(%) 

F330W 

EVPA(°) 

<*o-uv 

1 

3.1 ±0.3 

-79.9 ±3.12 

4.5 ±0.7 

-82.5 ±4.7 

1.50 ± 0.06 

40.2 ± 4.0 

-65.8 ±3.0 

34. 0 a ± 3.5 

-67.5 ±3.1 

0.07 ±0.04 

2 

3.7 ± 0.4 

-79.6 ±3.17 

4.5 ±0.7 

-82.5 ± 4.7 

1.41 ± 0.04 

39.9 ± 4.0 

-65.3 ±3.0 

34.0 ±3.5 

-67.5 ±3-1 

0.12 ±0.04 

3 

5.5 ± 0.6 

-70.3 ± 3.07 

.... 


1.23 ±0.08 

39.9 ± 4.0 

-51.3 ±3.0 



0.58 ±0.03 

4 

1.4 ±0.2 

-10.4 ±3.37 

13.4 ±1.5 

73.7 ±3.4 

1.54 ± 0.06 

23.4 ±2.3 

-58.2 ± 3.0 

24.8 ±2.6 

-64.7 ±3.1 

0.32 ± 0.03 

5 

2.0 ± 0.3 

-55.4 ±3.69 



1.58 ± 0.04 

27.7 ±2.8 

-56.9 ±3.0 



0.57 ±0.02 

6 

2.2 ± 0.3 

-98.9 ±3.76 



1.42 ± 0.04 

36.7 ±3.7 

-58.4 ±3.0 



0.66 ± 0.02 

7 

6.1 ±0.6 

-53.2 ±3.1 

.... 


1.41 ± 0.04 

42.5 ±4.3 

-61.5 ±3.0 



0.57 ±0.02 

8 

4.1 ±0.4 

-78.9 ±3.05 

10.0 ± 1.3 

-86.8 ± 3.8 

1.42 ± 0.04 

34.5 ± 3.5 

-64.7 ±3.0 

38.1 ±3.9 

-62.7 ± 3.0 

0.37 ±0.02 

9 

10.7 ± 1.1 

-36.0 ±3.04 



1.37 ±0.04 

32.4 ±3.2 

-63.8 ± 3.0 



0.35 ± 0.02 

10 

8.6 ±0.9 

-62.0 ±3.09 



1.48 ±0.05 

36.4 ± 3.7 

-60.6 ±3.0 



0.68 ±0.02 

11 

6.6 ±0.7 

-69.1 ±3.03 

10.1 ± 1.4 

-77.2 ± 4.0 

1.45 ± 0.05 

32.6 ±3.3 

-62.2 ± 3.0 

29.3 ±3.0 

-62.0 ± 3.0 

0.45 ±0.03 

12 

4.4 ±0.5 

-90.0 ±3.26 

.... 


1.44 ±0.05 

28.1 ±2.8 

-62.4 ±3.0 



0.67 ± 0.02 

13 

2.5 ± 0.3 

-69.4 ±3.24 

.... 


1.39 ±0.07 

26.3 ±2.6 

-59.5 ± 3.0 



0.55 ±0.03. 

14 

2.6 ±0.3 

-71.4 ±3.7 



1.50 ± 0.05 

24.1 ±2.4 

-59.0 ± 3.0 



0.63 ±0.03 

15 

1.0 ±0.1 

-1.2 ±3.91 

3.0 ± 1.0 

-25.7 ±9.8 

1.63 ±0.05 

20.7 ±2.1 

-62.7 ±3.0 

24.4 ± 2.7 

-50.2 ± 3.3 

0.47 ± 0.03 

16 

6.6 ±0.7 

-39.1 ±3.01 

2.0 ± 0.7 

-57 ± 10 

1.50 ±0.04 

34.9 = 3.5 

-66.1 ±3.0 

34.8 ± 3.6 

-64.2 ±3.1 

0.48 ±0.03 

17 

10.2 ± 1.0 

-46.4 ± 3.04 



1.60 ±0.04 

20.3 ±2.1 

-58.7 ± 3.0 



0.73 ± 0.03 

18 

6.7 ± 3.0 

-130.1 ± 1.30 

.... 


..... 

22.2 ±3.0 

-76.5 ± 0.7 

.... 




a Observations taken Dec. 10, 2002. 

two concentric apertures for the core (one with a radius 
of 11 pixels, the other with a radius of five pixels) and 
one centered on HST-1 with a radius of 11 pixels. The 
smaller of the two core apertures excludes a knot that 
emerged from the core in later images, while the larger 
core aperture includes it. While the two lightcurves show 
the same behavior, we use only the smaller aperture in 
this paper. For the lone WFPC2 dataset (epoch 18) we 
used apertures of 4 pixels radius for both the core and 
HST-1, and also excluded in each aperture an annular 
region between 4-6 pixels from the other. This alternate 
procedure was made necessary by the much larger pixels 
of the WFPC2/WFC chips. For the F606W images, af- 
ter galaxy subtraction rectangular regions located quasi- 
randomly throughout the galaxy were used to verify the 


flatness of the remaining background. A similar strat- 
egy was used for the F330W images, where the galaxy 
contributon was minimal. Aperture correction was done 
to account for the wide wings of the HST point spread 
function (PSF). As the nucleus and HST-1 are either un- 
resolved or nearly so, we modeled them as PSFs, using 
TINYTIM (Maybhate et al. 2010; Bohlin 2007a). The 
aperture corrections were generally 10-20% in flux, simi- 
lar to those found by Sirianni et al. (2005); however, our 
method allowed us to account for the fact that the posi- 
tion of these components varied widely on the chip from 
observation to observation. Count rates were converted 
to flux using SYNPHOT, recomputing PHOTFLAM val- 
ues by assuming a spectral index a = 0.7 (F u cx i/ -a ), ap- 
propriate for most of the M87 jet (Perlman et al. 2001). 
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Typical errors for this procedure are < 5%. Finally, we 
also accounted for Galactic extinction using data from 
NED, which gives Ab = 0.096, as well as standard ex- 
tinction curves. The resulting fluxes in Stokes I are listed 
in Table 1, while the fractional polarizations and EVPA 
are given in Table 2. 

High-quality UV photometry of these two components 
was published recently by Madrid (2009). We have uti- 
lized those measurements, albeit with some modificar 
tions. Rather than correcting all fluxes to a single wave- 
length as in the Madrid paper, we have utilized the un- 
corrected fluxes in both F220W and F250W. This mini- 
mizes the number of assumptions, and also allows us to 
make use of all the data in the 6 epochs where observa- 
tions were taken in both bands, increasing the accuracy 
of the measured spectral indices in those epochs and also 
allowing us to check for significant emission in the Mg 
IT A 2798 line, which falls near the center of the F250W 
bandpass, but is outside the F220W bandpass (no ev- 
idence of this line emission was found). The reader is 
referred to that paper for reduction steps. We list the 
resulting fluxes in Table 1. In Table 2 we list the spec- 
tral indices oluv-o- 

Errors were propagated in both datasets using stan- 
dard techniques. The propagated errors include Poisson 
errors, an additive noise term (the RMS background cal- 
culated post galaxy subtraction) and the read-out and 
discretization noises, as well as an additional 1% term to 
account for uncertainties in SYNPHOT models (Bohlin 
et al. 2007b). For this analysis, we have ignored errors in 
the F606W galaxy model. The resulting uncertainties in 
the Stokes Q and U images are approximately Gaussian 
in nature, with their values being approximately equal 
to the sum in quadrature of the individual polarizer im- 
age errors. Hence, Gaussian error propagation for P is 
appropriate for our purposes. 

4. RESULTS 

Figures 1 and 2 show our results for the total flux, 
opticai-UV spectral index, fractional polarization and 
EVPA variations of both the nucleus and HST-1. The 
total fluxes shown are in the F606W band, while all 
other panels use all available data. The total flux varia- 
tions mirror those previously shown for the near-UV by 
Madrid (2009), showing the very large flare in HST-1 
as well as two smaller flares in the nucleus during the 
same time. Because the large majority of these obser- 
vations occurred during 2004-2006, when the knot was 
already very bright, we do not see the full dynamic range 
of the variability exhibited by knot HST-1 during its 
flare - nearly a factor 150 at 2500 A, where its flux in- 
creased from 6 //Jy in 1999 to a peak of 854 fi Jy in early 
2005. We used our data to measure the distance between 
the two features, which in all epochs is consistent with 
0.885 = 0.010 arcsec, with no evidence of motion. This 
sets an upper limit of 1 .56 pc (2cr) on positional change 
of the flux maximum of HST-1 during the roughly 5-year 
timespan covered by these data. 

Comparing Figures 1 and 2 we see that the nucleus 
and HST-1 display very different characteristics in po- 
larization and spectral index variability. Both display 
large changes in polarization characteristics during the 
six-year timespan of our observations. The nucleus is 
seen to range between 1-13% polarization during this 


1 800 i- 
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FlG. 1. — Variations in Total flux (i.e. , Stokes I), fractional po- 
larization, spectral index and EVPA, axe plotted for the nucleus of 
M87. A dashed line in the EVPA panel reflects the jet PA. The 
F606W observations are plotted as squares, while in the second 
and fourth panel the F330W polarimetry is plotted as diamonds. 
See §3.2 and 4 for discussion. 
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Fig. 2. — Variations in Total flux, fractional polarization, spectral 
index and EVPA, are plotted for knot HST-1. A dashed line in 
the EVPA panel reflects the jet PA. The F606W observations are 
plotted as squares, while in the second and fourth panel the F330W 
polarimetry is plotted as diamonds. See §3.2 and 4 for discussion. 


time, with the EVPA changing by as much as 90 de- 
grees. HST-1 is much more highly polarized, with its po- 
larization ranging from 20-45%, but much less variability 
(marginally significant at most) in the EVPA, which in 
our data ranges from roughly —50° to —75°, with a typ- 
ical value of ~ -62°, about 7-8 degrees away from the 
nominal jet direction of —69.5° but rather closer (4 de- 
grees) to the average PA of the radius vector for HST-1 
during VLB A observations (Paper IV) of -66°. The dif- 
ference between the two angles is only 1 a for any one 
point, but in 17 of 18 epochs the EVPA is significantly 
displaced towards the north from the radius vector, mak- 
ing the difference significant at roughly the 3a level (see 
also §5.1). 

For the most part the polarization characteristics in 
F330W track the ones seen in F606W. This is true for 
all epochs for knot HST-1. Thus that knot displays no 
evidence for frequency-dependent polarization behavior. 
For the nucleus, however, one epoch shows significant (> 
2a) differences between polarization characteristics in the 
two bands, namely 28 Nov 2004. As can be seen in Figure 
1, this difference is highly significant both in P (8 a) as 
well as EVPA (10 a), and we have eliminated all possible 
sources of instrumental error in the F330W and F606W 
data for this epoch. We note that this epoch is very near 
the peak of a flare in the nucleus; however, this flare does 
not exhibit a spectrum that is significantly different than 
surrounding epochs. Unfortunately the next epoch does 
not have UV polarimetry so we are unable to comment 
further on whether any frequency-dependent polarization 
pattern developed during this flare. 

Another perspective can be gained by looking at rela- 
tionship between flux and fractional polarization. This 
is done in Figure 3, using the F606W data only. In 
knot HST-1 (bottom panel), this reveals a strong cor- 
relation between total flux and fractional polarization, 
particularly between 2004 November and 2006 December 
(epochs 4-17), when the flux variability was dominated 
by the main part of the flare and the monitoring was 
most frequent. The other four points represent times 
where the flux variability was either dominated by or 
had a significant contribution from smaller scale events. 
Using only epochs 4-17, a Spearman’s rho-test indicates 
p = 0.842 and P = 1.6 x 10~ 4 . Inspection of the F330W 
data shows they also follow this correlation. During this 
time the EVPA is nearly constant (Figure 1), although 
there is possible evidence for a quasi-sinusoidal modu- 
lation. This gives the clear impression of the variability 
being dominated by a single component with a highly or- 
dered magnetic field. We discuss the implications further 
in §5.2. 

A pattern is much more difficult to pick out for the 
nucleus. We have therefore added red arrows to Fig- 
ure 3 to guide the eye during the part of the campaign 
where there is evidence of a pattern. During epochs 6- 
15, there was a coherent increase and then decrease in 
the fractional polarization of the nucleus, at a time when 
the flux variations where small (Figure 1). This trans- 
lates into a near-vertical “loop” in the polarization-flux 
plane (Figure 3, top panel), with the initial increase in 
polarization being correlated with the decrease in flux 
from epoch 5 J s maximum, and then a nearly monotonic 
decrease being observed over epochs 9-15 while the flux 
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varied by only 25%. This is discussed further in §5.3. 
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FIG. 3.— Graphs of fractional polarization versus flux for the 
core (top) and HST-1 (bottom). The observation epochs have been 
labelled sequentially. Notice that for HST-1, epochs 4-18 display a 
very strong correlation between the flux and polarization, whereas 
the behavior for the nucleus is very different, with a ’loop’ seen 
between epochs 9-15, but otherwise no organized pattern. 

Table 2 and Figures 1-2 also describe the evolution of 
the UV-optical spectral index auv-o • As can be seen, 
&uv-o behaves differently for the nucleus than it does 
for HST-1. For the core, we do not see significant vari- 
ability in the optical spectral index. By contrast, the 
spectral index is strongly variable for HST-1. For most of 
the time, the variability of the spectral index appears un- 
correlated with flux, but during the brightest part of the 
flare (epochs 4-9, denoted on Figure 4 by red arrows), we 
see that the spectral index is larger (i.e., steeper) when 
the flux is higher. As can be seen in Figure 4, during 
this time HST-1 describes a definite ’loop’ in the flux- 
spectral index plane. This pattern has been called “hard 
lags” in the blazar and optical variability literature (see 
e.g., Zhang 2002, Zhang et al. 2002; Fiorucci, Ciprini & 
Tosti 2004; the term “counterclockwise looping” is also in 
use, but note that those papers use the opposite sign con- 
vention for ol than we do). Looping in the other direction 
is seen for epochs 13-17 (denoted on Figure 4 by green 
arrows), e.g., ’’soft lagging” (Zhang et al. 2002), with 


oscillations in the plane in between (denoted on Figure 
4 by blue arrows). Both “looping” patterns are known 
to arise for particular relations between the acceleration 
and cooling timescales controlling the spectral evolution 
of the radiating particles (in the framework of simplified 
models for such an evolution) and are discussed further 
in §5.2. 



Fig. 4. — The evolution of Q[/v — o versus flux for knot HST-1. 
Each epoch has been labelled sequentially, as in Figure 3. Notice 
the strong ’looping’ behavior during the maximum of the flare in 
HST-1, as well as after (green and red arrows). See §§4, 5.1 for 
discussion. 

5. DISCUSSION 

It is interesting to explore further the reasons behind 
the highly disparate behaviors of the nucleus and HST-1, 
as seen in our data as well as other bands. X-ray vari- 
ability of the nucleus and HST-1 was studied in Papers 
I, III and V. Variability has also been seen at gamma- 
ray energies, where M87 has been detected in both TeV 
(Aharonian et al. 2003) and GeV (Abdo et al. 2009) 
energies (although, n.b., in the gamma- rays the M87 nu- 
cleus and jet cannot be resolved from one another). Dur- 
ing the peak of the flare of HST-1 in March-May 2005, 
TeV variability on timescales of a few days (Aharonian 
et al. 2006) was seen, leading Stawarz et al. (2006) to 
associate the enhanced TeV emission seen in 2005 with 
the flare of HST-1. However, the origin of other features 
in the TeV lightcurve is unclear. A major work compil- 
ing all multi-wavelength and TeV gamma-ray variability 
data for M87 was recently completed by Abramowski et 
al. (2011). That work concluded that while it remains 
plausible that both the unresolved nucleus and HST-1 
contribute to the TeV emission observed from M87 sys- 
tem during the quiescence epochs and also the 2005 flare, 
during the 2008 and 2010 epochs of the enhaced 7-ray 
activity of the source the nucleus is more likely to have 
contributed the majority of the TeV flux. 

In the radio, both Paper IV as well as Chen et al. 
(2011) studied the variability of knot HST-1, with Paper 
IV using VLB A data, while Chen et al. (2011) used data 
from the VLA, both from roughly 2003-2007. The lat- 
ter work has angular resolution comparable to HST and 
includes an analysis of polarization data, a.nd finds a vari- 
able P , EVPA and rotation measure, as well as a radio 
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spectrum that softened when the knot was brighter. It is 
difficult to compare their data to ours in detail because 
Chen et al. (2011) used a subset of the available data, 
including only 4 datasets between 2004-2006 (and none 
in 2005), when HST-1 was most active. However their 
findings, while consistent with the idea of non-cospatial 
radio and optically emitting particle populations (Perl- 
man et al. 1999), are difficult to reconcile with the much 
longer radiative lifetimes of radio synchrotron emitting 
particles, as well as the similar radio and UV lightcurves. 

To further explore the physical implications of our re- 
sults, it is necessary to discuss the physics of shocks and 
other disturbances in jets, which may explain the be- 
havior we see. As will be seen, we do this because we 
believe that both behaviors may be explained by such 
disturbances. Following this, we will then present phys- 
ical interpretations of the behaviors of knot HST-1 and 
the nucleus. 

5.1. Shocks, Helical Distortions and Polarization 
Variability 

The behaviors we see in both the nucleus and HST-1 
are direct reflections of the physics in the emitting region. 
Since the data we have analyzed in this paper includes 
fluxes, optical- UV spectra and also polarimetry, we have 
information both on the interplay between particle accel- 
eration and cooling, as well as the magnetic field struc- 
ture that was either associated with this behavior or pro- 
duced it. In order to motivate the discussion herein, it is 
useful to summarize the commonalities in the behavior of 
both regions and then discuss why those commonalities 
argue for an origin in shocks and waves. Both HST-1 and 
the nucleus exhibit coherent patterns in the (7, P) plane. 
In the case of HST-1 the pattern is simple: polarization 
is correlated strongly with intensity, while at the same 
time the EVPA remains essentially constant very close 
to the PA of the jet. By comparison, in the nucleus we 
see a ’loop’ in the (7, P) plane, with somewhat more com- 
plicated EVPA behavior, featuring wild swings of up to 
100°, albeit around a dominant orientation that is once 
again near the PA of the jet’s radial motion vector. The 
fact that both HST-1 and the nucleus display coherent 
variability patterns in (I,P) while maintaining a single, 
dominant orientation of EVPA indicates that in both re- 
gions, the details of the local magnetic field structure are 
tightlv related to the efficiency of the particle accelera- 
tion. That component must have a reasonably well or- 
dered magnetic field structure, particularly in the case of 
HST-1 because of its very high polarization (see §5.2 for 
further elaboration on this issue). These behaviors are 
all consistent with having been produced through shocks 
and/or waves, although the details of the physics may be 
different in the two regions. 

What types of shocks may be consistent with these be- 
haviors? Perhaps the simplest type of shock to discuss 
is localized, planar, and oriented along or near the jet 
perpendicular. This type of model, often known as a 
“Laing sheet” due to the fact that shocks of this type 
characteristically compress an initially random magnetic 
field into a thin sheet with magnetic field along the sheet, 
has been investigated extensively in the literature, par- 
ticularly in Laing (1980), Hughes, Aller & Ailer (1985), 
and Kollgaard et al. (1990). In this model, the proper- 
ties of a relativistic shock can be completely determined 


by a few factors. Primary among these is its compres- 
sion ratio, k = TdPd/{V u Pu) (Laing 1980; Wardle et al. 
1994), where P Ui d is the bulk Lorentz factor and (3 u ,d is 
the bulk speed with the subscripts referring to upstream 
and downstream quantities respectively. One can also 
think of the compression ratio in terms of the compres- 
sion of a unit length. Ignoring the plasma magnetiza- 
tion, the shock jump conditions relate the upstream and 
downstream speeds via 0 u ftd =1/3 (Landau & Lifshitz 
1987) for a plasma with a fully relativistic equation of 
state. The other factors that describe the shock are the 
spectral index a, Doppler factor 6 and viewing angle 0 O &. 
More specifically, in this model, Kollgaard et al. (1990) 
found that the degree of polarization is given by (their 
equation (2); note that our convention for the spectral 
index a is the opposite of the one in Kollgaard et al.) 

3 + 3a <5 2 (1 - k 2 ) sin 2 9 ob 
5 4- 3a 2 — <5 2 (1 - k 2 ) sin 2 9 ob 

This model is obviously dependent on the geometry cho- 
sen for the jet magnetic field and for the shock - in par- 
ticular ignoring any helical component to the field and 
a strong, perpendicular shock. Such a model can eas- 
ily produce a correlation between flux and polarization, 
along with a roughly constant EVPA parallel to the jet 
direction, as seen in knot HST-1. We will discuss in §5.2 
the application of this model to HST-1. 

If, however, a different type of disturbance is envi- 
sioned, quite different polarization behavior can be pro- 
duced. The first possibility we will consider is a time 
varying upstream speed, fi u . This is supported by the 
variable super luminal speeds observed in AGN jets(e.g. 
Lister et al. 2009). In this second model, the varying 0 U 
is related to the downstream flow speed via fi u fid = 1/3 
which induces variability in the shock compression fac- 
tor, k. To induce a variation in the shock compression 
factor, the stochastic variation in fi u is modeled as being 
sinusoidal: 

Pu{t) = Pu,o + Ap sin ( ujpt - <pp) (2) 

The parameters Ap,ujp, 4>p) are free parameters of 
the model, and t is time. If the jet contains a disordered 
field component B r , no ordered toroidal field, and an or- 
dered poloidal field B p , then the polarization depends 
on the magnetic field through the ratio £ = \B P /B r \. 
Since the process of transforming an initially ordered 
large-scale field to a tangled one is probably governed 
by the kink instability (Spruit, Foglizzo & Stehle 1997; 
Begelman 1998; Marscher 2009), it is natural to assume 
that the value of £ at the shock fluctuates in time. Sim- 
ilar to p u , we model the £ variation as 

€(t) = So + A ? sin (u^£ - 0 4 ). (3) 

The parameters (£o> u/^, fa) are free parameters of the 

model, and t is time. Equations (2) and (3) allow the cal- 
culation of the post shock intensity, 7, and polarization 
P (e.g. Kollgaard et al. 1990): 

/ oc 5 2 d +a k~ 2 B 2 T {2 + [3 kY - (1 - fc 2 )] sin 2 9' ob ) (4) 
p „ 3 + 3a 8 2 [(1 - k 2 ) - 3fc 2 £ 2 ] sin 2 9 ob 

5 — 3a 2 — (5 2 (1 — fc 2 ) sin 2 9 ob + 3<5 2 £ 2 fc 2 sin 2 9 ob 

(5) 
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Note that equation (4) is an approximation as it results 
from integrating over the line of sight with a set to unity 
for convenience (Wardle et al. 1994), and that equation 
(5) as expected reduces to equation (1) for £ = 0. The 
downs'.ream Doppler factor, 6 = (r </ - Td0dcos9 ob )~ 1 , 
raised to the power 2 + a consistent with a steady jet 
(Lind & Blandford 1985), is included in the expression 
for intensity since the flow speed varies in the emission 
region which is downstream of the shock. In our conven- 
tion, the EVPA is parallel to the jet axis for P > 0, and 
is perpendicular to it for P < 0. The factor B 2 in equa- 
tion (4) is held as constant in our model. Otherwise, if 
used as another time varying parameter, it would merely 
modify the amplitude of the intensity fluctuations. 

Another plausible configuration is one where the jet 
contains a significant helical component to its mag- 
netic field. Such a scenario might be particularly op- 
erative in the nucleus, since AGN jets are generally 
thought to be launched and collimated within the cen- 
tral parsec by magnetic fields that dominate other 
sources of pressure. Its applicability to regions further 
from the nucleus is less certain because it is not known 
how far downstream from the launching site the jet re- 
tains a magnetically dominant ordered helical field since 
jets are thought to be unstable to the m = 1 kink 
mode (Narayan et al. 2009; Marscher 2009; Spruit 
2010). Despite this theoretical uncertainty, some obser- 
vational evidence suggests that jets contain helical mag- 
netic fields on parsec scales: parsec-scale bulk acceler- 
ation (Vlahakis & Konigl 2004), Faraday rotation gra- 
dients (Asada et al. 2002; Gabuzda, Murray & Cronin 
2004; Zavala & Taylor 2005; Kharb et al. 2009), 
and observed asymmetries in the transverse pro- 
files of polarization, brightness and spectral index 
(Clausen-Brown, Lyutikov & Kharb 2011). 

Let us assume that a component in the jet contains a 
large-scale helical magnetic field that is variable due to 
the growth of the m = 1 kink mode. The jet’s magnetic 
symmetry axis (parallel to the B z direction) will then 
deform into a large scale helix that is carried with the 
jet’s velocity field (Mizuno et al. 2011), which we as- 
sume to be uniform in this work. If the helical deforma- 
tion passes through a standing shock, as shown in Figure 
5, then the magnetic structure of the post-shock region 
will change in time, thereby producing fluctuations in the 
post-shock synchrotron emission. Let us further assume 
that the ratio of B'JB^ varies in time in the post-shock 
flow. Kink mode simulations have found that, to avoid 
total jet disruption, B^/B' dynamically relaxes to ~ 1 
(Nakamura et al. 2007), while, in competition with this 
process, jet conical expansion always increases the ratio. 
Therefore, this competition will produce fluctuations in 
B'JB'z at the standing shock. 

Within this scenario multiple types of disturbances can 
be envisioned. Herein we consider two. Firstly, we may 
introduce sinusoidal variations both in the jet frame mag- 
netic pitch angle ip' = tan _1 (B£/B'), and in the angle 
between the magnetic symmetry axis and the jet frame 
line of sight, 9' , as shown in Figure 5 : 

ib'(t) = A ^ sin (uty£ - <£,/ ) + ip' 0 (6) 

0'{t) = A' 0 sin (cj<?£ - 4> e ) + 0' ob (7) 


> 

amission region 



Fig. 5. — Schematic of the helical distortion model illustrated 
with toroidal loops of magnetic field. As shown by the velocity 
vectors (labeled by v), the velocity field is uniformly directed de- 
spite the kinked jet. As the kinked region of the jet propagates 
through the standing shock’s emission region (gray region), the 
orientation of the magnetic field with respect to the line of sight 
changes with time. This change in orientation is parametrized as 
6'{t) in equation (7). The effects of shock compression, the poloidal 
magnetic field and the change in the magnetic pitch angle, ip 1 (t) , 
are not shown. 

where t is time and (A' e , are 

parameters of the model. The jet frame viewing an- 
gle, 0' ob , is actually set by the relation to the observer 
frame viewing angle by sin0^ 6 = <5sin0o6. However, as 
the Doppler factor for the inner jet is unconstrained, d' ob 
is treated as a free parameter as long as the required 
Doppler factor is within reasonable bounds. It should be 
noted that when the magnetic field passes through the 
standing shock, the field components lying in the shock 
plane will be amplified by shock compression. However, 
this only modifies the form of ip'(t ) and 9'{t)\ it does not 
prevent quasi-periodic variations from occurring. 

To calculate the total intensity and fractional polarize 
tion from the standing shock, the emission region elec- 
tron distribution function is assumed to be a power law, 
dn = K e E~ p dE, where the spectral index is related to 
the electron distribution function by p = 2a 4- 1. We 
assume the jet is unresolved and the emission is mostly 
concentrated in a cylindrical shell centered on the local 
symmetry axis so that the total intensity and fractional 
polarization are (Lyutikov, Pariev & Gabuzda 2005): 

I & K (cos 2 ip' + cos 2 9' - 3(cos 9 1 cos ip') 2 -I- 1) (8) 

p ~ 3 + 3a -2(1 + 3 cos 2iP') sin 2 9' 

5 + 3a 5 — cos 29' — cos 2ip' — 3 cos 29' cos 2 ip' ’ 

( 9 ) 
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where O' and ip' are the variables defined in equa- 
tions (7) and (6) respectively, and if is an arbi- 
trary constant that depends on emission region details 
such as the beaming factor, emission region size, rel- 
ativistic particle density, and magnetic field strength 
(Lyutikov, Pariev & Gabuzda 2005). The sign of P in- 
dicates whether the EVPA is parallel (P > 0) or perpen- 
dicular (P < 0) to the jet frame local magnetic symmetry 
axis in the standing shock. At any time t 1 , the values of O' 
and ip' represent the particular orientation of the kinked 
magnetic field and the value of the magnetic pitch angle 
in the standing shock respectively. (See the appendix for 
a derivation of equations 8 and 9.) 

5.2. Polarization and Spectral Behavior of HST-1: 
Interpretation 

We believe that the most consistent explanation for the 
variability observed in HST-1 is that the flare occurred 
in a shock within the jet, with the maximum polariza- 
tion coming at the time of maximum compression and 
also maximum optical flux. The high polarization (P at 
maximum in excess of 40%) and alignment of the EVPA 
with the jet axis in HST-1, rules out the conical recon- 
finement hydrodynamic shock model of Nawaiejko (2009; 
see also Bromberg & Levinson 2009). Assuming that 
the jet magnetic field is weak and tangled, such mod- 
els cannot produce a polarization higher than ~ 30%. 
Furthermore, such a model would predict a substantially 
different orientation of EVPA for small and intermediate 
jet inclinations. The EVPA we observe in HST-1 is very 
nearly perpendicular to the jet, and nearly constant, so a 
more consistent explanation for the data is that the non- 
thermal activity is restricted to a localized, perpendicular 
(possibly stationary) strong shock within the interior of 
the flew, as envisioned in §5.1 (see eq. (1)), which is 
also more in line with the observed radio morphology. 
The feature producing the flare may then be associated 
with the Mach disk produced around the nozzle (con- 
verging point) of the reconfinement shock (Stawarz et 
al. 2006, Paper IV, Gracia et al. 2009). In the latter 
model (Gracia et al. 2009) the synchrotron emissivity 
is proportional to the comoving frame electron density 
and magnetic field strength and configuration, modulo 
£( 2+a ). rather than being dependent on hydrodynamic 
energy dissipation as assumed in' Nawaiejko (2009). 

Using the doubling/halving timescales calculated for 
the optical/UV emission in Madrid (2009) as well as 
those in Paper II for the X-ray emission, we can estimate 
a size of < 0.5b light-years (where <5 is the Doppler fac- 
tor) for the size of the flaring region, with the irregularity 
of the available measurements perhaps arising because of 
the complexity of the compression mechanism (we do not 
believe it arises because we did not sample adequately, 
given the smoothness of the observed lightcurve). This 
is consistent with the fact that HST-1 is not resolved 
by HST, which sets a hard upper limit of 2 pc on the 
radius of the optically emitting components (Paper II, 
Madrid 2009). Interestingly, however, as already noted 
in Paper II and Madrid (2009), the X-ray and optical/UV 
doubling/halving timescales do not lead to significantly 
different constraints for the region size. 

If we then apply equation (1) to knot HST-1, we can 
reproduce the type of behavior seen in the (/, P) plane 
(Figure 3). But equally interestingly, we can also con- 


strain the kinematics in the local flow under certain as- 
sumptions. To illustrate, we choose k = 0.25, appropri- 
ate for a strong relativistic shock (Meisenheimer et al. 
1989), spectral index 0.45 < a < 0.55, corresponding to 
the mean observed during both flare “loops” (Figures 2, 
4), and 0.35 < P < 0.45, consistent with the peak value 
attained at both flux maxima. We then allow the view- 
ing angle 0 o b and Doppler factor S to vary. We constrain 
the solutions so that the apparent speed 0 aPl) falls within 
the range 4.0 < /3 app < 4.5, the values reported for VLBI 
components in Paper IV. The result is shown in Figure 
6. As can be seen, the permitted values of 0 o i range be- 
tween 11 — 18°, and 6 can range between 2.5 and 5.5. 
The Lorentz factor, T is more stable, however, with per- 
mitted values ranging between 4.1 and 4.8. These values 
are in agreement with the upper limit of S = 8 calculated 
by Waters & Zepf (2005) based on other considerations, 
and are also consistent with the requirement of 0 o b < 19° 
from the somewhat faster superluminal speeds seen in 
HST monitoring during the 1990s (Biretta et al. 1999). 
It should be mentioned that with these values of a , P 
and k there are no allowed solutions with T < 2, and 
moreover, nearly all the allowed solutions with T < 3 
require values of 0 O b > 20°, which is not allowed. 

Interestingly however, we see no evidence for motion of 
the flaring region of HST-1 in our data, as our astromet- 
ric results (all epochs having identical distances between 
the nucleus and HST-1 to a tolerance of ±0 01") translate 
to an upper limit of 1.02 c (2cr) on motion of the flaring 
region itself. While this seems to conflict with the VLBA 
measurements of Paper IV, as well as our calculation of 
5 from the shock model, this should not be too concern- 
ing. It is entirely possible that the flaring region itself 
represents a standing shock, not flowing along with the 
plasma. Under such a scenario we could not expect the 
observed speed to accord with the local jet Lorentz fac- 
tor. Indeed, the VLBA maps (Paper IV) show a station- 
ary upstream end to the HST-1 region. For most of the 
time period covered by Paper IV, the observed speeds are 
consistent with our limit, as the flaring region (compo- 
nent HST-lc in their nomenclature), while downstream 
from the stationary upstream end (HST-ld), has a speed 
of 1.14 ± 0.14 c , measured relative to HST-ld. Begin- 
ning around epoch 2006.0, however, HST-lc splits into 
two components, with the faster, downstream one (ac- 
counting for the majority of the radio flux) accelerating 
to 4.3 ±0.7 c relative to HST-ld. While this latter speed 
is on the surface highly inconsistent with our astrometric 
results, it is important to realize that our data are not 
very sensitive to this time period, as it contains only 5 of 
our 17 ACS epochs (the lower-resolution WFPC2 obser- 
vations of epoch 18 are much less useful for astrometry). 
Using only our ACS data during 2006, the limit on mo- 
tion of the flaring region would be much less restrictive, 
i.e., ~ 5.5 c at 2 o. Use of later epoch data would improve 
this result; however, there was a gap in UV monitoring 
of the M87 jet between late 2006 and mid-2009, a time 
interval that featured a further factor ~ 5 decrease in 
the X-ray flux of HST-1 (e.g., Paper V and later data), 
which could make it very difficult to identify the flaring 
components. The necessary data do, however, exist in 
the VLBA archive, and it would be highly interesting to 
track the motion of the flaring region in both bands. 

If indeed the giant HST-1 flare was a result of en- 
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FlG. 6. — Allowed values of beaming parameters for our shock 
model of knot HST-1. At top, we show the allowed range for the 
Lorentz factor V and viewing angle 0 o b, while at bottom, we show 
the corresponding range for the Doppler factor 5 plotted against 
viewing angle 0 O & . In both panels, light gray colors refer to spectral 
indices a; = 0.45, gray colors refer to a = 0.5, and black colors refer 
to a = 0.55 See §§5.1, 5.2 for discussion. 


hanced particle acceleration occurring within a compress- 
ing shock, then it would be logical to ascribe the cool- 
ing during the main flare to relaxation of the compres- 
sion that occurred within the shock. In that case, the 
energy dependence - or lack thereof - of the variabil- 
ity timescales (in optical and other bands) may be set 
by the dynamical timescale of the compression, as origi- 
nally noted in Paper II and explored further in Paper V 
and Madrid (2009). For example, if the compression and 
subsequent expansion were adiabatic, one would expect 
to see frequency-independent variability behavior unless 
there was an intrinsic, pre-existing break in the spec- 
trum. This would be modified where the particle cooling 
and acceleration timescales (£ occ and t coo i respectively) 
are equal to or less than the dynamical timescales. As al- 
ready noted, we observe hard lags in the optical-UV dur- 
ing the brightest part of the flare (epochs 5-9, Figure 4). 
To explain such hard lags, t acc must be similar to t coo i, 
whereas when soft lagging was observed (epochs 13-17), 
the opposite relationship would hold, i.e., t acc < tcooi- 


Interestingly, as shown in Paper V, in both the UV and 
X-rays the derivative dl /dt changed sign between 2005.4- 
2005.5, i.e., epochs 9 and 10. If indeed this was related 
to the relationship between the acceleration and cool- 
ing timescales, then the oscillations seen in epochs 10-13 
- which occurred during a time when the flux was de- 
creasing monotonically - become important. The mul- 
tiwavelength spectral characteristics of HST-1 discussed 
above could be possibly explained in the framework of 
the scenario in which the primary loss mechanism in the 
optical-UV is Comptonization of external radiation and 
the optical-UV emitting electrons are near the transition 
between the Thomson and Klein-Nishina regimes (see the 
discussion in §6, below). 

As discussed in §4, the mean EVPA in HST-1 is some- 
what different from the PA of the jet. Thus while a shock 
morphology is likely for HST-1, the polarization data 
hints at a more complex morphology than indicated by 
the simple, unresolved appearance shown by the images. 
Two interpretations are possible. The first possibility is 
that HST-1 ’s optical polarized flux comes from two or 
more regions. While this might seem the simplest inter- 
pretation, it is difficult to reconcile this with the strong 
correlation between flux and polarization, combined with 
the constant EVPA. In addition, the unresolved nature of 
the flaring region in our data and that of Madrid (2009) , 
constrains the maximum separation of these components 
to ~ 2 parsecs, which translates to a constraint on the 
cooling timescale that, as already discussed, is consis- 
tent with the spectral evolution we see. Alternatively, 
it could indicate either a twist in the field within the 
shocked region and/or an oblique shock, as suggested for 
the knot A region by Bicknefl & Begelman (1996). In 
the latter case, one might expect to see evidence of a 
slight local deviation in the flow direction, and indeed, 
the VLBA components seen in Paper IV do have a sig- 
nificant range of radius vector PAs as measured from 
the standing feature at the upstream end of the HST-1 
complex. Furthermore, there is weak evidence of small 
changes in EVPA over time, with the EVPA near the 
two flux maxima (epochs 8, 9 and 16) being within 2-3 
degrees of the PA of the radial vector from the nucleus 
to HST-1 from the VLBA data, while in the other 15 
epochs the EVPA is more closely aligned with the mo- 
tion vector of the faintest feature of the HST-1 complex 
as seen on the VLBA maps, namely component “a”, at 
the downstream end. Note however that if the shock lies 
well within the interior of the jet these deviations might 
not be indicative of an overall deflection of the flow (as in 
knot A). On a somewhat related note, Nakamura et al. 
(2010) recently proposed a model for the overall struc- 
ture of the M87 jet, in which the main shocks in the M87 
jet have a helical magnetic field structure, which could be 
produced via two methods. Bicknell & Begelman (1996) 
suggested that the knots in the M87 jet are generated by 
the helical modes, which can produce lateral oscillation of 
the entire jet but has trouble producing filaments within 
it. Hardee & Eilek (2011), by contrast, postulate that the 
elliptical Kelvin-Helmholtz modes dominate, and model 
the inner M87 jet as having twisted, high-pressure fil- 
aments generated by elliptical Kelvin-Helmholtz insta- 
bilities which eventually disrupt the jet beyond knot A. 
Assuming that the magnetic field within the components 
is perpendicular to the local bulk velocity vector (rather 



Polarization and Spectral Variability in M87 


11 


than the mean jet axis vector), such a model can produce 
the slight offset between EVPA and jet PA seen in nearly 
all of the observations of HST-1. 

5.3. Polarization Behavior of the Nucleus: 
Interpretation 

Because the polarization and flux variability of the nu- 
cleus forms a coherent pattern-a counter-clockwise loop 
in I - P space, over the course of a year-the avail- 
able data challenge scenarios where the variability of the 
nucleus is dominated by multiple independent compo- 
nents. The most natural interpretation of the nuclear 
polarization variability is that within the small aperture 
chosen (which represented approximately 2 HST resolu- 
tion elements) the structure of the M87 jet is relatively 
simple, perhaps having a single region of the jet dom- 
inating the flux and polarization variability at any one 
time. This scenario is consistent with either models of 
a magnetically dominated jet in which large-scale mag- 
netohydrodynamic instabilities play an important role 
(e.g. Giannios & Spruit 2006) or with models of a non- 
stationary flow through a standing shock. Within such a 
model there are a number of ways to explain the behavior 
seen in the nucleus. In §5.1 we have illustrated two pos- 
sible scenarios that may occur - namely, either (1) there 
is a standing shock through which a time varying flow 
propagates, or (2) the jet contains a large-scale helical 
field subject to the kink mode in a region of a standing 
shock as shown in Figure 5. Both of these possibilities 
(see the discussions surrounding equations (2)-(5) and 
(6)-(9) respectively), are motivated by (i) the variability 
of jet flow speeds as seen in superluminal studies of jets 
and (ii) the prominence of the kink mode in theoretical 
analyses of magnetically dominated jets respectively. 

The two scenarios considered here each produce differ- 
ent tracks in the (/, P) plane. In Figure 7, we show 
an example track for each, assuming a viewing angle 
6 0 h = 15°, consistent with the observation of superlu- 
minal motions at speeds as high as 6c in HST-1 (Biretta 
et al. 1999, Paper IV; see also §5.2). As shown, both sce- 
narios successfully reproduce two general characteristics 
of the looping behavior we see in the (/, P) plane be- 
tween epochs 6 to 15 in Figure 3, namely: (i) successive 
points are located close together in (/, P) space, and (ii) 
these points form a coherent pattern (a loop) over the 
course of a ~ year as shown in figure 7. Nonetheless, 
both models suffer from shortcomings. Constraining the 
shock compression model to produce EVPAs parallel to 
the jet direction severely restricts the maximum polariza- 
tion it can achieve. Thus, the model’s maximum polar- 
ization reaches only ~ 6%, while the observations achieve 
a maximum of ~ 11%. In the helical distortion model, 
maintaining a low polarization and producing an EVPA 
parallel to the jet implies restricting the variability of the 
magnetic pitch angle to be ip' = 1/2 cos -1 (-1/3) -I- e, 
where 0 < e <C 1, which is an arbitrary constraints on 
the oscillation of ip 1 . More generally, both models are 
limited by their extreme simplicity (e.g. neither have 
sheared velocity fields), high number of free parameters, 
and non- uniqueness. 

Of the two toy models discussed herein, we would ar- 
gue that the (/, P) behavior we observe (Figure 3) is 
closer to that produced by the helical distortion model, 
although again, neither model reproduces the observed 



Intensity (arbitrary units) 

Fig. 7. — Illustrated here is a theoretical counter-clockwise I - 
P loop for the helical distortion model (solid) and the variable 
shock compression model (dotted). The values used in the shock 
compression model are: {&v,O>A 0 ,U 0 ,<pp) = (0.78, 0.03, 1, — tt/ 2) 
and ({o,^,«{,^) = (0.78, —0.5, 1,0). For the helical distortion 
model, {A' g ,u}e,(j)e,Q' ob ) — (1.5, l,0,7r/2) and (A(,. , , tp' 0 ) = 

(0.05,2,0,1.03) and <5 = sin0' 6 /sin 15° « 3.86. For both models, 
9 a h = 15° and a = 1.5 were used. 

behavior perfectly. The shock compression model re- 
quires a roughly linear increase of polarization with flux 
in its ’positive’ stage, followed by a ’plateau’ with po- 
larization nearly constant at its maximum while the flux 
decreases. Neither of these is observed in the nucleus. In- 
stead we see very steep increases in polarization during 
epochs 6-9 (from 2% up to 12%), while the flux changes 
are much smaller. The pattern does not exactly repro- 
duce what is seen in the helical distortion model either, 
however. Local maxima in flux are seen at epochs 7 and 
12. Of these, the former corresponds to a local maxi- 
mum in P while the latter does not; indeed, the global 
maximum in P comes at epoch 9, when the flux is de- 
clining. In fact, one could actually argue that a ’plateau’ 
is seen between epochs 9 and 10, as the two polarization 
measurement are statistically indistinguishable from one 
another while the fluxes differ by 25%. It is possible that 
both types of variations are seen, somewhat out of phase 
with one another, but while this is plausible such a model 
would produce competing variations that are difficult to 
model without introducing ad hoc assumptions. 

Neither of these models is fully successful in reproduc- 
ing the EVPA fluctuations we observe (Figure 1). The 
helical distortion model predicts that the EVPA time 
evolution should correlate with the polarization and in- 
tensity. We expect this correlation because the helical 
distortion flowing through the standing shock causes the 
orientation of the magnetic field symmetry axis in the 
post shock flow to fluctuate as illustrated in figure 5; in 
turn, the post-shock flow EVPA, intensity, and polarize 
tion depend on the fluctuating field geometry. Unlike the 
helical distortion model, the varying shock compression 
model maintains cylindrical symmetry even as other pa- 
rameters vary with time. Thus, the shock model predicts 
the EVPA should remain stationary or undergo rapid 90° 
flips. It should be noted, however, that the dominant ob- 
served EVPA of the nucleus is close to parallel to the jet, 
an orientation which both toy models reproduce for the 
parameters used in figure 7. 

The lack of variability we observe in ao-uv for the 
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nucleus ostensibly fits better within the helical distor- 
tion than with the shock compression model. In the 
helical distortion model, the shock strength and parti- 
cle distribution parameters, K e and p (= 2a + 1), are 
held constant and only the magnetic geometry fluctuates, 
hence the lack of spectral index variation in the model. 
However, more realistically, in particle acceleration mod- 
els such as diffusive shock acceleration (or first-order 
Fermi acceleration), the particle acceleration efficiency 
depends on the angle between the large-scale magnetic 
field and the shock plane (Sironi & Spitkovsky 2011). In 
the case of the shock compression model, the increased 
shock compression ratio leads to enhanced particle accel- 
eration, as envisioned in §5.2 for HST-1, the expectation 
is a strong correlation between flux and auv-o, either 
hard or soft lagging depending on the relationship be- 
tween the acceleration and cooling timescales. 

6. SUMMARY 

This work has revealed the value of polarization and 
spectral information in monitoring campaigns, as well as 
in understanding the physics of jet regions. We have 
shown that knot HST-1 can be best understood as a 
shocked region, where the flaring upstream end did not 
move significantly between 2002-2007, and displays char- 
acteristics that are consistent with a classic perpendicu- 
lar shock. Under such a model we find that the Lorentz 
factor in the jet at HST-1 can be constrained to a fairly 
narrow range, 4.1 < Y < 4.8, and the viewing angle can 
also be constrained to 10° < 0 O & < 18°. The optically 
emitting part of the jet within HST-1 may well contain 
multiple components (as the radio emitting part does, see 
Paper IV), given the fine points of the EVPA structure 
(§5.2); however the HST data does not have the angular 
resolution to resolve this region. The spectral behavior of 
HST-1 was consistent with nearly equal particle acceler- 
ation and cooling timescales in the optical- UV; however, 
the fact that the X-ray emission is best understood as 
synchrotron radiation complicates things. Two interpre- 
tations are possible for this duality. The first of these is 
that the X-ray emission may come from a small part of 
the jet in HST-1 (as suggested for downstream regions 
of the jet by Perlman & Wilson 2005), which may either 
be partially co-spatial with the optical emission region or 
distinct from it. This interpretation is disfavored, how- 
ever, because of the overall similar appearance of the 
UV and X-ray lightcurves (Paper V). A more likely in- 
terpretation is that, while the likely emission mechanism 
in both bands remains synchrotron radiation, the dom- 
inant energy loss mechanism is inverse-Comptonization 
of external radiation, with the X-ray emitting electrons 
being in the Klein-Nishina regime, while the optical/UV 


emitting electrons are either in the Thomson regime or 
around the Thomson/Klein-Nishina transition. Under 
such a scenario the cooling timescale of lower-energy 
optical-UV photons may be comparable to or shorter 
than the cooling timescale of X-ray emitting electrons. 
The flat optical spectra observed in HST-1 (a £ 0.5) are 
in fact very consistent with such an idea (Moderski et 
al. 2005). In such a scenario the inverse-Compton emis- 
sion would dominate at higher energies and in fact would 
be energetically dominant when integrated over the en- 
tire electromagnetic spectrum. This would be consistent 
with the correlation of the 2005 TeV flare with the HST-1 
flare seen at lower energies (Abramowski et al. 2011). 

The variability behavior seen in the nucleus, however, 
can be best understood as either a helical distortion to 
a steadv jet, where the distortion would arise from kink 
mode instabilities, or fluctuations in the jet speed that 
produce corresponding fluctuations in the strength of 
shocks within the nuclear jet. Both of these toy mod- 
els can produce the fluctuations in field components that 
can result in the ’looping’ polarization behavior we see. 
We believe the most likely mechanism is a current-driven 
instability combined with a fluctuating, helical magnetic 
pitch angle. Jet precession can be ruled out as there 
is no evidence of jet wobbles on larger resolved scales. 
This leaves large-scale instabilities as a possible cause of 
the wobbling axes. Kelvin- Helmholtz instabilities could 
cause the wobbling; however these are important in hy- 
drodynamic jets and as discussed in Section 5.3 we be- 
lieve that in the nuclear regions the M87 jet is more likely 
to be magnetically dominated, with a helical magnetic 
field structure. The kinked jet scenario naturally leads 
to magnetic pitch angle fluctuations. Future monitor- 
ing of the variability behavior of spatially resolved jets 
should include both multi-band imaging and polarimetry 
on the same timescale in order to maximize the physical 
information that can be gained from the campaign. 
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APPENDIX 


HELICAL FIELD INTENSITY AND POLARIZATION 

Here we derive the expressions for synchrotron intensity (eqn. 8) and polarization (eqn. 9) of an unresolved jet with 
an emitting cylindrical shell containing a helical field, a configuration first considered in the nonrelativistic case by 
(Laing 1981). The notation and discussion here follows that of Lyutikov, Pariev & Gabuzda (2005), who explicitly 
derive the polarization expression in their equation (21). The synchrotron intensity and polarization can be expressed 
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as integrals over the entire cylindrical shell: 

I = K |B , sinx'r fl <ty 
Jo 

3 + 3a / 0 27r \B' sin x! |“ +1 cos 2 xd<p 


P = 


5 + 3a f 2n \B f sm X '\ a+l &t> 


(Al) 

(A2) 


where <p is the azimuthal cylindrical coordinate, K is a constant, x! is the angle between the jet frame line of sight 
and the jet frame magnetic field, and x is the angle between the EVPA and the jet’s projected onto the sky, measured 
clockwise. To obtain analytical results, we assume a = 1 in carrying out the above integration. 

Assume a Cartesian coordinate system centered on the jet with the bulk flow directed along the 2 -axis and with the 
observer in the y = 0 plane. Therefore, the jet velocity is 


and the photon propagation vector is 


/3 = /?(0,0,l) 


h = (sin B ob , 0, cos 9 ob ) 


(A3) 


(A4) 


The jet frame shell magnetic field unit vector has a pitch angle of tan = B'^j B' z and can be expressed in cylindrical 
coordinates as 


B' = (— sin ip' sin <£, sin ip' cos (p , cos ip') . 


(A5) 


Thus sin xJ —n' • B' can now be expressed as: 

sin 2 x! = cos 2 ip' sin 2 B' ob + ^ sin 2 8' ob sin 2 ip' sin <p + (cos 2 B' oh + cos 2 <p sin 2 9 ob ) sin 2 ip' (A6) 

Integrating this expression over cp for a = 1 produces the intensity 

I = K (cos 2 ip' + cos 2 0' ob — 3(cos B' ob cos ip') 2 + l) ( A7) 

To find the polarization, all that remains to be calculated is cos 2^. First, note that the polarization vector of a 
synchrotron electromagnetic wave in the jet frame is e! = n' x B' . x may now be written as 

cos x' = e' ■ (n' x gj , (A8) 

where ft' x £ is the jet direction projected onto the sky, and £= (0, 1,0) in Cartesian coordinates. Evaluating cosx' 
leads to 


cos 2x' = 


cos 2 <p sin 2 ip' - (cos ip' sin 9' ob cos 0' ob sin <p sin ip') 


i\ 2 


Therefore, according to equation (A2), 

P = 


1 — (cos ip' cos 0' ob — sin B' ob sin (p sin ip')‘ 

3/2 (1 + 3 cos 2V’ 7 ) sin 2 6’ oh 
5 - cos 2 6' ob - cos 2 ip' — 3 cos 9' ob cos 2 ip' ' 


(A9) 


(A10) 


Note that equations (A6) and (A9) are equations (20) of Lyutikov, Pariev & Gabuzda (2005). Although this result 
strictly holds for a = 1, it can be extended to other values close to a = 1 by the following analytic approximation 


/ 3 + 3a \ 2 (1 + 3 cos 2V> 7 ) sin 2 9' ob 

\ 5 + 3a / 5 - cos 2 6' ob - cos 2ip' - 3 cos 9' ob cos 2 ip' ' 


(All) 
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